Artificial spin-orbit coupling in ultra-cold Fermi superfluids 



Kangjun Seo, Li Han and C. A. R. Sa de Melo 
School of Physics, Georgia Institute of Technology, Atlanta, Georgia 30332, 

(Dated: October 31, 2011) 



USA 



The control and understanding of interactions in many particle systems has been a major chal- 
lenge in contemporary science, from atomic to condensed matter and astrophysics. One of the 
most intriguing types of interactions is the so-called spin-orbit coupling - the coupling between the 
spin (rotation) of a particle and its momentum (orbital motion), which is omnipresent both in the 
macroscopic and microscopic world. In astrophysics, the spin-orbit coupling is responsible for the 
synchronization of the rotation (spinning) of the Moon and its orbit around Earth, such that we 
can only see one face of our natural satellite. In atomic physics, the spin-orbit coupling of electrons 
orbiting around the nucleus gives rise to the atom's fine structure (small shifts in its energy levels). 
In condensed matter physics, spin-orbit effects are responsible for exotic electronic phenomena in 
semiconductors (topological insulators) and in superconductors without inversion symmetry. Al- 
though spin-orbit coupling is ubiquitous in nature, it was not possible to control it in any area of 
physics, until it was demonstrated in a breakthrough experiment [l[| that the spin of an atom could 
be coupled to its center-of-mass motion by dressing two atomic spin states with a pair of laser beams. 
This unprecedented engineered spin-orbit couplingwas produced in ultra-cold bosonic atoms, but 
can also be created for ultra-cold fermionic atoms [J-yl- In anticipation of experiments, we develop 
a theory for interacting fermions in the presence of spin-orbit coupling and Zeeman fields, and show 
that many new superfluids phases, which are topological in nature, emerge. Depending on values 
of spin-orbit coupling, Zeeman fields, and interactions, initially gapped s-wave superfluids acquire 
p-wave, d-wave, /-wave and higher angular momentum components, which produce zeros in the 
excitation spectrum, rendering the superfluid gapless. Several multi-critical points, which separate 
topological superfluid phases from normal or non-uniform, are accessible depending on spin-orbit 
coupling, Zeeman fields or interactions, setting the stage for the study of tunable topological super- 
fluids. 

PACS numbers: 03.75.Ss, 67.85.Lm, 67.85.-d 



The effects of spin-orbit coupling in few body systems 
like the Earth-Moon complex in astrophysics or the elec- 
tron spin and its orbital motion around the nucleus in 
isolated atoms of atomic physics are reasonably well un- 
derstood due to the simplificity of these systems. How- 
ever, in the setting of many identical particles, spin-orbit 
effects have revealed quite interesting surprises recently 
running from topological insulators in semiconductors [J| 
to exotic superconductivity Q and non-equillibrium ef- 
fects depending on the precise form of the spin-orbit 
coupling. In atomic physics the coupling arises from the 
interaction of the magnetic moment of the electron and 
a magnetic field, present in the frame of electron, due 
to the electric field of the nucleus. Similarly in con- 
densed matter physics, the coupling arises from the mag- 
netic moment m of electrons, which move in the back- 
ground of ions. In the electron's reference frame, these 
ions are responsible for a magnetic field B, which de- 
pends on the electron's momentum k and couple to elec- 
tron's spin. The resulting spin-orbit coupling has the 
form Hso ~ — m • B = —^jhj{k)aj, where CTj rep- 
resents the Pauli matrices and /ij(k) describes the j-th 
component [j = x, y, z) of the effective magnetic field 
vector h. For some materials h can take the Dressel- 
haus 01 form hij(k) = voikyX + fc^-y), the Rashba [1| 
form h/{(k) = vii{—kyX + kxy), or more generally a lin- 
ear combination of the two h^(k) = h£i(k) + hj^(k). In 
all these situations the type of spin-orbit coupling can 



not be changed arbitrarily and the magnitude can not 
be tuned from weak to strong, making the experimental 
control of spin-orbit effects very difficult. 

Recently, however, it has been demonstrated experi- 
mentally that spin-orbit coupling can be engineered in 
a ultra-cold gas of bosonic atoms in their Bose-Einstein 
condensate phase [l[ , when a pair of Raman lasers creates 
a coupling between two internal spin states of the atoms 
and its center-of-mass motion (momentum). Thus far, 
the type of spin-orbit field that has been created in the 
laboratory [l| has the equal-Rashba-Dresselhaus (ERD) 
form h^(k) = hERoi^) = vk^y, where vr = vd = v/2. 
Other forms of spin-orbit fields require additional lasers 
and create further experimental difficulties In ultra- 
cold bosons the momentum-dependent ERD coupling has 
been created in conjunction with uniform Zeeman terms, 
which arc independent of momentum, along the z axis 
(controlled by the Raman coupling ^r), and along the 
y-axis (controled by the detuning 5). The simultaneous 
presence of hz, hy and IierdO^) leads to the Zeeman- 
spin-orbit (ZSO) Hamiltonian 
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for an atom with center-of-mass momentum k and spin 



basis It), U). The fields 



~^r/2, hy = -5/2 and 



^ERD ~ vkxH can be controlled independently, and thus 
can be used as tunable parameters to explore the avail- 
able phase space and to investigate phase transitions, as 
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achieved in the experiment involving a bosonic isotope of 
Rubidium («^Rb). Ahhough current experiments have 
focused on Bose atoms, there is no fundamental reason 
that inipeeds the realization of a similar set up for Fermi 
atoms [ll-Q designed to study fermionic superfluidity Q . 
Considering possible experiments with fermionic atoms 
such as ^Li, "^^K, we discuss in this letter phase diagrams, 
topological phase transitions, spectroscopic and thermo- 
dynamic properties at zero and finite temperatures dur- 
ing the evolution from BCS to EEC superfluidity in the 
presence of controllable Zeeman and spin-orbit fields in 
three dimensions. 

To investigate artificial spin-orbit and Zccman fields in 
ultra-cold Fermi superfluids, we start from the Hamilto- 
nian density 



7{(r) =Ho(r)+H/(r), 
where the single-particle term is simply 



(1) 
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Here, Ka = — V^/ (2m) — /io, is the kinetic energy in refer- 
ence to the chemical potential , hj (r) is the combined 
effective field including Zccman and spin-orbit compo- 
nents along the j-dircction (j = x,y,z), and iptii^) arc 
creation operators for fermions with spin a at position 
r. Notice that we allow the chemical potential to be 
different from /ij^, such that the number of fermions N-^ 
with spin ^ may be different from the number of fermions 
with spin ^. The interaction term is 



Hiir) 



gVt(r)V^I(r)V;;(r)i/'t(r), 



(3) 



where g represents a contact interaction that can be 
expressed in terms of the scattering length via the 
Lippman-Schwinger relation V/g ~ -~Vm/{4:Tras) + 
l/(2ek)- The introduction of the average pairing field 
A(r) = g{tl)i{r)ip^{r)) « Aq and its spatio-temporal fiuc- 
tuation 77(r, t) produce a complete theory for superfluid- 
ity in this system. 

From now on, we focus on the experimental case where 
a) the Raman detuning is zero {S = 0) indicating that 
there is no component of the Zccman field along the y di- 
rection; b) the Raman coupling flu is non-zero meaning 
that a Zeeman component along the z direction exists, 
that is, hz = —Vlfi/2; and c) the spin-orbit field has com- 
ponents hyCk) and /ia;(k) along the y and x directions. To 
start our discussion, we neglect fluctuations, and trans- 
form Ho(r) into momentum space as i/o(k). Using the 
basis V'j(k)|0) = |k t), V'j(k)|0) = |k 1), where |0) is 
the vacuum state, the Fourier-transformed Hamiltonian 
-ffo(k) becomes the matrix 

Ho(k) = K+(k)l + K.az - h^cr^ - hy{k)ay - ha;(k)(Tx, 

Such matrix can be diagonalized in the helicity basis 
$[(k)|0) EEE |k $]^(k)|0) = |k 4), where the spins 



and -IJ. are aligned or antialigned with respect to the 
effective magnetic field hoff(k) = h||(k) + hj^(k). Here, 
K+Ck) = (A'-f + Ki)/2 = Ek — M+7 is a measure of the 
average kinetic energy = k'^/2m in relation to the 
average chemical potential fi+ = {fj,^ + ni)/2. While 
hx(k) = /ia,(k)x -I- hy (k)y is the spin-orbit field and 
h||(k) = {hz — K^)z is the effective Zeeman flcld, with 
K- = {Kf - Ki)/2 = where = (^^ - /^^)/2 is 
the internal Zeeman field due to initial population im- 
balance, and is the external Zeeman field. When 
there is no population imbalance the internal Zccman 
field is fi- = 0, and we have only hz. In general, 
the eigenvalues of the Hamiltonian matrix Ho(k) are 
e^(k) - A + (k)-|h ff(k)| ande^W ^ A-+ (k) + |heff(k)|, 
where |hcH(k)| = a/ -I- h^)^ + |/i_L(k)p is the magni- 
tude of the effective magnetic field, with the transverse 
component being expressed in terms of the complex func- 
tion h±{'k) = hx{k.) + ihy{\i). In the limit where the in- 
ternal H- and external hz Zeeman fields vanish and the 
spin-orbit field is null {h± = 0), the energies of the helic- 
ity bands are identical 5ff(k) = £ji.(k) producing no effect 
in the original energy dispersions [ifll • 

When interactions are added to the problem, pair- 
ing can occur within the same helicity band (intra- 
helicity pairing) or between two different helicity bands 
(inter- helicity pairing). This leads to a tensor order 
parameter for superfluidity that has four components 
Ajf^(k) = — AT(k)e^"^, corresponding to the helicity 
projection A = -1-1; A^4(k) = -A5(k), and A||.^(k) = 
A5(k), corresponding to hehcity projection A = 0; and 
A4|4(k) = — Ar(k)e*''', corresponding to helicity pro- 
jection A = —1. The phase (p(k) is defined from 
the amplitude-phase representation of the complex spin- 
orbit field h^{k) = \h^{k)\e"^^^\ while the amplitude 
AT(k) = Ao|/i±(k)|/|hcff (k)| for helicities A = ±1 are di- 
rectly proportional to the scalar order parameter Ag and 
to the relative magnitude of the spin-orbit field |/ij_(k)| 
with respect to the magnitude of the effective magnetic 
field |hoff(k)|. Additionally, At has the simple physi- 
cal interpretation of being the triplet component of the 
order parameter in the helicity basis, which is induced 
by the presence of non-zero spin-orbit field h± , but van- 
ishes when h± = 0. Analogously the amplitude A5(k) = 
Ao/i|| (k)/|hoff (k)| for helicity A = are directly propor- 
tional to the scalar order parameter Aq and to the rela- 
tive magnitude of the total Zeeman field (k) = + hz 
with respect to the magnitude of the effective magnetic 
field |heff(k)|. Additionally, As has the simple physical 
interpretation of being the singlet component of the order 
parameter in the helicity basis. It is interesting to note 
the relation |AT(k)|2 + |As(k)|2 = |Aop, which, for fixed 
|Ao|, shows that as |A5'(k)| increases, |AT(k)| decreases 
and vice-versa. Such relation indicates that the singlet 
and triplet channels arc not separable in the presence of 
spin-orbit coupling. Furthermore, the order parameter in 
the triplet sector A|f^(k) and AD.jj.(k) contains not only 
p-wave, but also /-wave and even higher odd angular mo- 
mentum contributions, as long as the total Zeeman field 
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li-+hz is non-zero. Similarly, the order parameter in the 
singlet sector A^jj.(k) and A4^(k) contains not only only 
s-wave, but also d-wave and even higher even angular 
momentum contributions, as long as the total Zeeman 
field /i_ + hz is non-zero. Higher angular momentum 
pairing in the helicity basis, occurs because the original 
local (zero-ranged) interaction in the original (f, i) spin 
basis is transformed into a finite-ranged interaction in 
the helicity basis (ff, JJ.). In the limiting case of zero total 
Zeeman field fi_ + hz = 0, the singlet component vanishes 
(A5(k) = 0), while the triplet component becomes inde- 
pendent of momentum (AT(k) = Aq), leading to order 
parameter Aj).^(k) = — ft,^(k), and Ajui.(k) = —h±{]i.) 
which contains only p-wave contributions [ll| , since the 
components of h± (k) depend linearly on momentum k. 

The eigenvalues Ej{\i) of the Hamiltonian including 
the order parameter contribution emerge from the diago- 
nalization of a 4 x 4 matrix (see supplementary material). 
The two eigenvalues for quasiparticles are 



E,{k) = y (6.- - + |As(k)Pj + |AT(k)P, (4) 

corresponding to the highest-energy quasiparticle band, 
and 



i?2(k) = ^(a_ + v'a^ + |As(k)p) +|AT(k)P, (5) 

corresponding to the lowest-energy quasiparticle band, 
while the eigenvalues for quasiholes are i?3(k) ~ — £'2(k) 
for highest-energy quasihole band and E4(k) = —Ei(k) 
for the lowest-energy quasihole band. The energy ^ 
[£,-fiOi-) — Cj|(k)] /2 is momentum-dependent, corresponds 
to the average energy difference between the helicity 
bands and can be written as = — |heff(k)|, while 
the energy = [^f|-(k) -I- ^4j,(k)] /2 is also momentum 
dependent, corresponds to the averaged energy sum of 
the helicity bands and can be written as ^/j^ — _fC+(k) — 
Ek - M+- 

There are a few important points to notice about the 
excitation spectrum of this system. First, notice that 
£'i(k) > E2(k) > 0. Second, that the eigenergies are 
symmetric about zero, such that we can regard quasi- 
holes (negative energy solutions) as anti-quasiparticles. 
Third, that only £'2(k) can have zeros (nodal regions) 
corresponding to the locus in momentum space satisfying 

the following conditions: a) ~ ^xj^h+ + |As(k)p, 

which corresponds physically to the equality between the 
effective magnetic field energy |heff(k)| and the excita- 
tion energy for the singlet component ^J^h^ + |As(k)p; 

and b) |AT(k)| = 0, corresponding to zeros of the triplet 
component of the order parameter in momentum space. 

Since -B2(k) < £'i(k), and only i?2(k) can have ze- 
ros, the low energy physics is dominated by this ein- 
genvalue. In the case of equal Rashba-Dresselhaus 
(ERD) where h±(k) = zeros of i52(k) can oc- 

cur when kx =0, leading to the following cases: (a) 



two possible lines (rings) of nodes at {ky + kl)/{2m) = 
fi+ + \/ {fi- + hzY — |Aop for the outer ring, and {ky + 
k^)/ (2m) = /i+ — ^/(JI~\~Jlz)^~-JAo\^ for the inner ring, 



|Ao| 



> 0; (b) doubly-degenerate 
/i+ for > 0, doubly- 



when (/z_ + hz)^ 
line of nodes {ky + k1)/{2m 
degenerate point nodes for /i+ = 0, or no-line of nodes 
for ^1+ < 0, when (^_ + hz)"^ - |Aop = 0; (c) no 
line of nodes when (/i_ -I- /i^)^ — |AoP < 0. In ad- 
dition, case (a) can be refined into cases (a2), (al) 
and (aO). In case (a2), two rings indeed exist pro- 
vided that /i+ > y^/tZ+7i7)2^^|A^. However, the 
inner ring disappears when fi^ ~ ^ (/i_ + /i^)^ — |AoP. 
In case (al), there is only one ring when j^+j < 
yj {\i- -I- ft,^)^ — |AoP, In case (aO), the outer ring dis- 
appears at [ij^ — —yj (/i_ + hz)^ — |Aop, and for /x+ < 
— \/{ji- + hzY — |Aop no rings exist. 

We choose our momentum, energy and velocity scales 
through the Fermi momentum kp^ defined from the to- 



tal density of fermions 



71^ 



+ n; = fcj,^/(37r2). This 



choice leads to the Fermi energy ep^ = kp^/2m and to 
the Fermi velocity vp^ = kp^/m, as energy and veloc- 
ity scales respectively. In Fig. 1, we show the phase 
diagram of Zeeman field hz/ep^ versus chemical poten- 
tial ^+/eF+ describing possible superfiuid phases accord- 
ing to their quasiparticle excitation spectrum. We la- 
bel the uniform superfiuid phases with zero, one or two 
rings of nodes as US-0, US-1, and US-2, respectively. 
Non-uniform (NU) phases also emerge in regions where 
uniform phases are thermodynamically unstable. The 
US-2/US-1 phase boundary is determined by the condi- 
tion i^L+ = v/(Ai- +hzy - |Ao|2, when + hz\ > \ Aa\; 
the US-O/US-2 boundary is determined by the Clogston- 
like condition |(/i_ + hz)\ = |Ao| when > 0, where 
the gapped US-0 phase disappears leading to the gap- 
less US-2 phase; and the US-O/US-1 phase boundary 
is determined by /i+ = — ■\/ (/i_ -I- /iz)^ — jAgp, when 
+ hz\ > |Ao|. Furthermore, with the US-0 bound- 
aries, a crossover line between an indirectly gapped and 
a directly gapped US-0 phase occurs at = 0. Lastly, 
some important multi-critical points arise at the inter- 
sections of phase boundaries. First the point ^+ = and 
|(^_ + hz)\ = |Ao| corresponds to a tri-critical point for 
phases US-0, US-1, and US-2. Second, the point |Ao| = 
and = \{fj.- + hz)\ corresponds to a tri-critical point 
for phases N, US-1 and US-2. In the limit where both 
/i_ and hz vanish no phase transitions take place and the 
problem is reduced to a crossover [T2| - [l^ . 

In the US-1 and US-2 phases near the zeros of i?2(k), 
quasiparticles have linear dispersion and behave as Dirac 
fermions. Such change in nodal structures is associated 
with bulk topological phase transitions of the Lifshitz 
class as noted for p-wave 15 1 and d-wave [3, [13 superfiu- 
ids. Such Lifshitz topological phase transitions are possi- 
ble here because the spin-orbit coupling field induces the 
triplet component of the order parameter AT(k). The 
loss of nodal regions correspond to annihilation of Dirac 
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FIG. 1: Phase diagram of Zeeman field h^/epj^ versus cliem- 
ical potential ^j^/ep^ for a) v/vp^ = and b) v/vp^ = 0.28 
identifying uniform superfluid phases US-0 (gapped), US-1 
(gapless with one ring of nodes), and US-2 (gapless with 
two-rings of nodes). The NU region corresponds to unsta- 
ble uniform superfluids which may include phase separation 
and/or a modulated superfluid (supersolid). Solid lines rep- 
resent phase boundaries, while the dashed line represents the 
crossover from the direct-gap to the indirect-gap US-0 phase. 

quasiparticles with opposite momenta, which lead to the 
disappearance of rings. The transition from phase US- 
2 to indirect gapped US-0 occurs through the merger of 
the two-rings at the phase boundary followed by the im- 
mediate opening of the indirect gap at finite momentum. 
However, the transition from phase US-2 to US-1 corre- 
sponds to the disappearance of the inner ring through 
the origin of momenta, similarly the transition from US- 
1 to the directly gappped US-0 corresponds to the dis- 
appearance of the last ring also through the origin of 
momenta. In the case of Rashba-only coupling rings 
of nodes are absent and it is possible to have at most 
nodal points [l^, [l^. The last two phase transitions 
are special because the zero-momentum quasiparticles at 
these phase boundaries correspond to true Majorana zero 
energy modes if the phase ^'(k) of the spin-orbit field 
/ij_(k) = |/i_L(k)|e*^('') and the phase e{k) of the order 
parameter Ag = |Ao|e*^('^' have opposite phases at zero 
momentum: (p{0) — —0{O) [mod(27r)]. This can be seen 
from an analysis of the quasiparticle eigenfunction 

$2(k) ^ Ul(k)7/'kt + W2(k)V'ki + U3(k)V'lkt + "4(k)V'lki 

corresponding to the eigenvalue i?2(k). The emergence 
of zero-energy Majorana fermions requires the quasiparti- 
cle to be its own anti-quasiparticle: $2('^) ^ 'I'2(k). This 
can only happen at zero momentum k = 0, where the am- 
plitudes ui(0) — 1X3(0) and 7/2(0) = u^iO). Such require- 
ment leads to the conditions fi^ ~ (/x_ + h,)'^ + |Aop, 
and (p{0) — —0{O) [mod(27r)], showing that Majorana 
fermions can exist only at the US-O/US-1 and US-2/US- 
1 phase boundaries. It is important to emphasize that 
the Majorana fermions found here exist in the bulk, and 
thus their emergence or disappeareance affect bulk ther- 



modynamic properties, unlike Majorana fermions found 
at the edge (surfaces) of topological insulators and some 
topological superfluids. The common ground between 
bulk and surface Majorana fermions is that both exist at 
boundaries: the bulk Majorana zero-energy modes may 
exist at the phase boundaries between two topologically 
distinct superfluid phases, while surface Majorana zero- 
energy modes may exist at the spatial boundaries of a 
topologically non-trivial superfluid. 

It is evident that the transition between different su- 
perfluid phases occurs without a change in symmetry in 
the order parameter Aq, and thus violates the symmetry- 
based Landau classification of phase transitions. In the 
present case, the simultaneous existence of spin-orbit and 
Zeeman fields (internal or external) couple the singlet 
A5(k) and triplet AT(k) channels and all the super- 
fiuid phases US-0, US-1 and US-2 just have different 
weights from each order parameter component. How- 
ever a finer classification based on topological charges 
can be made via the construction of topological invari- 
ants. Since the superfluid phases US-0, US-1, US-2 are 
characterized by different excitation spectra correspond- 
ing to the eigenvalues of the Hamiltonian matrix includ- 
ing interactions H(k), we can use the resolvent matrix 
R(w, k) =_[— wl + H(k)] and the methods of algebraic 
topology [20[ to construct the topological invariant 

where Afc^^ = R9fc^,R^^. The topological invariant is 
£ = in the gapped US-0 phase, is £ = 1 in the gap- 
less US-1 phase and ^ = 2 in the gapless US-2 phase, 
showing that, for ERD spin-orbit coupling, £ counts the 
number of rings of zero-energy excitations in each super- 
fluid phase. The integral above has a hyper-surface mea- 
sure dSj and a domain 2? that encloses the region of zeros 
of a; = EjCk) = 0. Here ^, ly, A, 7 run from to 3, and 
has components ko = OJ, ki = kx, k2 = ky, and = k^. 
The topological invariant measures the flux of the four- 
dimensional vector i^T = e'^'^^TTr [Afe^,Afc,>feJ /247r2, 
through a hypercube including the singular region of the 
resolvent matrix R(a;, k), much in the same way that the 
flux of the electric field E in Gauss' law of classical elec- 
tromagnctism measures the electric charge q enclosed by 
a Gaussian surface: § dS ■ (eqE) = q. Thus, the topologi- 
cal invariant defined above defines the topological charge 
of fermionic excitations, in the same sense as Gauss' law 
for the electric fiux defines the electric charge. 

A full phase diagram can be constructed only upon 
verification of thermodynamic stability of all the pro- 
posed phases. For this purpose it becomes imper- 
ative to investigate the maximum entropy condition 
(see supplementary material). Independent of any 
microscopic approximations, the necessary and suffi- 
cient conditions for thermodynamic stability of a given 
phase are: positive isovolumetric heat capacity Cy = 
T {dS/dT)y ^^^-^ > 0; positive chemical susceptibility 
matrix ^0,^3 = {djjLa/ dNi3)rp y , i.c, eigenvalues of the 
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matrix [^] are both positive; and positive bulk mod- 
ulus B — or isothermal compressibility kt = 
— V~^ {dV/dP)j. ^j^^^ . Using these conditions, we con- 
struct the full phase diagram described in Fig.[l]for equal 
Rasha-Dresselhaus (ERD) spin-orbit coupling. The re- 
gions, where the uniform superfluid phases are unsta- 
ble are labeled by the abbreviation NU to indicate that 
non-uniform phases such as phase separation or modu- 
lated superfluid (supcrsolid) may emerge. In Fig. [2l we 
show the phase diagram of Zeeman field hz/ep^ versus in- 
teraction parameter l/{kp^as), for population balanced 
fermions, where the number of spin-up fermions iV-j- is 
equal to the number of spin-down fermions iVj,. 
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FIG. 2: Phase diagram of Zeeman field hz/ep^ versus interac- 
tion 1/ {kp+CLs) showing uniform superfluid phases US-0, US- 
1, and US-2, and non-uniform (NU) region for a) v/vf_^ ~ 0; 
b) v/vF+ = 0.14; c) v/vF+ = 0.28; d) v/vf^ = 0.56. 
Solid lines are phase boundaries, the dashed line indicates 
a crossover from the indirect- to direct-gapped US-0. 

Since these superfluid phases exhibit major changes in 
momentum-frequency space as evidenced by their single 
particle excitation spectrum, it is important to explore 
additional spectroscopic quantitities to characterize fur- 
ther the nature of these phases and the phase transitions 
between them. An important quantity is the 4x4 resol- 
vent matrix 

/ G(ia;,k) F(iw,k)\ 
^ ' yFT(iw,k) G(ia;,k) J ^ ' 

from where the spectral density ^o,(a;,k) = 
— (l/7r)ImGaQ(ia; — uj + iS,k) for spin a =t,i can 
be extracted. The spectral function ylct(a;,k) in the 
plane of momenta ky-k^ with fc^ = and frequency 
a; = reveals the existence of rings of zero-energy 
excitations in the US-1 and US-2 phases. The density of 
states = ^j^^Q(a;,k) for spin a as a function of 

frequency w is also an important spectroscopic quantity 



which is shown in Fig. [3] along with excitation spectra 
Ej{\C} for phases US-1 and US-2 at fixed ERD spin-orbit 
coupling v/vpj_ — 0.28. The parameters used for phase 
US-1 are hz/tp^^ = 0.5 and l/(fcF+as) = —0.4, while for 
phase US-2 they are hz/ep^ — 2.0 and l/{kp^as) — 1.0. 
Notice that, even though the excitation spectrum Ej{\^) 
is symmetric, the coherence factors appearing in the 
matrix G are not, such that the density of states 2?q(w) 
is not an even function of uj, and thus it is not particle- 
hole symmetric. The main feature of I?Q,(a;) at low 
frequencies is the linear behavior due to the existence 
of Dirac quasiparticlcs and quasiholes in the US-1 and 
US-2 phases, which are absent in the direct-gap and the 
indirect-gap US-0 phases. The peaks and structures in 
X'a(cj) mostly emerge due to the maxima and minima of 
Ej{k). Notice that for finite Zeeman field hz, the density 
of states 2?-|-(cj) ^ '^^{lxj) because the induced population 
imbalance P = {N^ — N{) / [N^^ + N]^) is non-zero. For 
the US-2 case shown in Fig. [3]d, the induced population 
imbalance P 1 since hz/ep^ is small, while for the 
US-1 case shown in Fig. |3j;, P ~ 1 as the spins are 
almost fully polarized since hz/tp^ is large. 



a b c 




/fe///CF+ D(a))£F^ Ikyl/kF^ 

d e f 




i l- I I.I -ri -U.BH I I I li -n.» h I I ■ I -H 

0.5 1 1.5 2 0.5 1 1.5 2 0.5 1 1.5 2 



IkJ/kF* D(co)£f^ Ikyl/kF^ 

FIG. 3: Energy spectrum and density of states in phase US-2 
are shown in a), b), c) for hz/ep^ = 0.5 and l/{kF_^_as) ~ —0.4 
and in phase US-1 are shown in d), e), f) for h^/ep^ = 2.0 
and l/(fcF+ffls) = 1.0. Energies Ej{kx, 0,0) versus \kx\ in a) 
and d); frequency uj versus density of states T>f{Lj) (dashed), 
'Di{iu) (dot-dashed), and their sum T>{i-o) (solid) in b) and e); 
energies Ej{0,ky,0) versus \ky\ in c) and f). 

In summary, we have discussed the effects of spin-orbit 
and Zeeman fields in ultra-cold Fermi superfluids, ob- 
tained the phase diagrams of Zeeman field versus inter- 
action parameter or versus chemical potential, and iden- 
tified several bulk topological phase transitions between 
gapped and gapless superfluids as well as a variety of 
multi-critical points. We have shown that the presence of 
simultaneous Zeeman and spin-orbit fields induces higher 
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angular momentum pairing, as manifested in the emer- 
gence of momentum dependence of the singlet and triplet 
components of order parameter expressed in the helicity 
basis. Finally, we have characterized topological phases 
and phase transitions between them through their exci- 
tation spectra (existence of Dirac quasiparticles or Majo- 
rana zero-energy modes), topological charges, and spec- 
troscopic and thermodynamic properties, such as density 



of states and isothermal compressibility. 
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Z = ^'D\^,i\^\ciL\){~S\^,il^\) with action 



tl - 



drdr 



^^i(r,T) — V'a(r,T)+H(r,T) 



where the Hamiltonian density is given in Eq. ([!]). 

Using the standard Hubbard-Stratanovich transfor- 
mation that introduces the pairing field A(r,T) = 
g{^Pl{r, T)ip^(r, r)) and integrating over the fermion vari- 
ables lead to the effective action 



SeS 



drdr 



|A(r,r) 



T , , M 
27 ^"'^"^Y 



K+S{r 



where K+ = (Kf + Ki)/2. The matrix M is 



M 



/ dr+ A't -h±^ 
-hi dr + A ^ 

At - K-^ 



\ -At 



(7) 







dr - ) 



where h±^ ~ — ihy corresponds to the transverse com- 
ponent of the spin-orbit field, to the parallel com- 
ponent with respect to the quantization axis z, — 



A'-f- — hz , and 



To make progress, we use the saddle point approxi- 
mation A(r,r) = Ao -|- 77(r,r), and write M = Mgp + 
Mf. The matrix Mgp is obtained via the saddle point 
A(r, r) Ao which takes M — )■ Mgp, and the fluctua- 
tion matrix Mf = M — Mgp depends only on 7y(r, t) and 
its Hermitian conjugate. Thus, we write the effective ac- 
tion as S'cff = Ssp + S{. The first term is 



ARTIFICIAL SPIN-ORBIT COUPLING IN 
ULTRA-COLD FERMI SUPERFLUIDS: 
(SUPPLEMENTARY MATERIAL) 



E 

]<.,iuJn -j 



In 



T 



k 



in momentum- frequency coordinates (k, iaj„); where 
ujn — {2n + 1)ttT. Here, Ej(k) are the eigenvalues of 



The method used to study the spin-orbit and Zeeman 
effects in ultra-cold Fermi superfluids is the functional 
integral method and its saddle-point approximation in 
conjunction with fluctuation effects. To describe the 
thermodynamic phases and the corresponding phase dia- 
gram in terms of the interactions, Zeeman and spin-orbit 
fields, we calculate partition function at temperature T 



/ ift(k) -^^±(k) 



H 



sp 



V 



-Ao 
-/jl(k) X^(k) ^Ao 

At -A>(-k) /zl(-k) 

-At /;a(-k) -i^i(-k) / 

(8) 

which describes the Hamiltonian of elemen- 
tary excitations in the four-dimensional basis 



7 



*t = {v;j(k),Vl(k),Vt(-k),V'i(-k)} 
ation action is 

■|7y(r,r)|2 T 



The fluctu- 



drdv 
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In dot (l + M-^Mf) 



The spin-orbit field is h^(k) = h/f(k) + h£)(k), 
where hij(k) = vn {—ky'k + A:3:y) is of Rashba-type and 
h£i(k) = vd {ky'k + kxf) is of Dresselhaus-type, has 

magnitude |/ij_(k)| = (un - vuf + (ud + vnf k^. 

For Rashba-only (RO) {vd = 0) and for equal Rashba- 
Dresselhaus (ERD) couplings {vu = vo ~ the mag- 
nitude of the transverse fields are |/i±(k)| = VR^k'^ + fc^ 

{vr > 0) and h±{'k) = v\kx\ (v > 0), respectively. 

The Hamiltonian in the helicity basis $ = U^', where 
U is the unitary matrix that diagonalizes the Hamilto- 
nian in the normal state, is 



H,p(k) 









a(k) 

A, 



Airt(k) 

AM(k) 



VA;^(k) A*^(k) 



AM(k) \ 

AM(k) 


-a(k) / 



The components of the order parameter in the helic- 
ity basis are given by A-ff^(k) = AT(k)e~"'''', and 
A^ll.(k) = — AT(k)e*''''' for the triplet channel and by 
A^4(k) = -As(k) and A4^(k) = As(k) for the sin- 
glet channel. The eigenvalues of Hsp(k) for quasiparti- 
clcs i?i(k), £'2(k) are hsted in Eqs. and ©, while 
the eigenvalues for quasiholes are E^Ck) = —E2(k), and 
Ei{k) = -Ei{k). 

The thermodynamic potential is = flgp + ^^f , where 
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T 



- ^ In {1 + exp [^E,{k)/T]} + Yl 



with K+ = [-ft^t(-k) + Ki{-k) /2 is the saddle 

point contribution and flf = — TlnZf, with Z{ = 
J 'D[fj,r]]exp[—S{{fj,r/)] is the fluctuation contribution. 
The order parameter is determined via the minimization 
of Qsp with respect to |Aop, leading to 



^=-lT.^F[E,{k)] 



aiAoP ■ 



(9) 



where np [Ej{k)] = l/(exp [£'j(k)/T] + 1) is the Fermi 
function for energy Ej(k). The contact interaction g is 
expressed in terms of the scattering parameter Os via the 
Lippman-Schwinger relation discussed in the main text. 

The total number of particles N+ = + is defined 
from the thermodynamic relation 7V_|_ = — [dU,/ dyb^)rp y , 
and can be written as 



'Sp T- iVf. 

The saddle point contribution is 



(10) 



sp 



+ / T,V 



1 - 



and the fluctuation contribution is iVf 
— {dilf/dfi+, )y y leading to 

N, = ^Jv% ^] exp [-S,{n, V)] 
with the partial derivative being 



T 
2V 



-Tr 



(1 + M-iMf) ' 



Knowledge of the thermodynamic potential ft, of the 
order parameter Eq. ([9|) and number Eq. (jlOp provides 
a complete theory for spectroscopic and thermodynamic 
properties of attractive ultra-cold fermions in the pres- 
ence of Zeeman and spin-orbit fields. Representative 
Saddle point solutions for chemical potential /x+ and or- 
der parameter amplitude |Ao| as a function of l/{kp^as) 
in the equal Rashba-Dresselhaus (ERD) case {v/vf^. = 
0.28) are shown in Fig. H for h^/eF+ = 0,0.5,1.0,2.0. 
These parameters are used to obtain the phase diagrams 
described in Figs. [1] and [2] in combination with an anal- 
ysis of the excitation spectrum Ej{k) given in Eqs. (|3]) 
and ([5|) and the thermodynamic stability conditions for 
all the uniform superfluid phases: directly or indirectly 
gapped superfiuid with zero nodal rings (US-0); gapless 
superfluid with one ring of nodes (US-1); and gapless 
superfluid with two rings of nodes (US-2). 




-10 12 
1/kF+as 



1 2 

1/kF+as 



FIG. 4: a) Chemical potential fi+/eF_^. and b) order pa- 
rameter amplitude |Ao|/eF_,, versus interaction parameter 
l/{kF^as) for spin-orbit parameter v/vp^ = 0.28 and val- 
ues of the Zeeman field h^/ep^ ~ (solid); h^/ep^ ~ 0.5 
(dashed); h^/ep^ — 1.0 (dotted); and h^/ep^ — 2.0 (dot- 
dashed). 

A thermodynamic stability analysis of all proposed 
phases can be performed by investigating the maximum 
entropy condition. The total change in entropy due to 
thermodynamic fluctuations, irrespective to any approx- 
imations imposed on the microscopic Hamiltonian, can 
be written as 



ASt, 



1 

2T 



(ATAS* - APAy + Afi^ANc,) , 



where the repeated a index indicates summation, and 
the condition A^tot l£ guarantees that the entropy is 
, maximum. Considering the entropy S" to be a function 
of temperature T, number of particles Na and volume 
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V, we can elliminate the fluctuations AS, AP, and A^q, 
in favor of fluctuations AT, Ay and AiV„, and show 
that the fluctuations AT are statistically independent of 
ANa and AV, while fluctuations ANa and AV are not. 
The first condition for thermodynamic stability leads 
to the requirement that the isovolumctric heat capacity 



Cv =T (dS/dT) 



> 0. Additional conditions are 



directly related to number AA^^ and volume AV fluctu- 
ations. They require the chemical susceptibility matrix 
= [dlJ'aldNp)j,y to be positive definite, i.e, that 
its eigenvalues are both positive. This is guaranteed by 
det[^] = Ctt^ii - ^tiCit > and ^tt > 0- The last con- 
dition for thermodynamic stability is that the bulk mod- 
ulus B = 1/ kt or the isothermal compressibility kt = 
— V^^ {dV/dP)rp J , are positive. Since the number 
ANa and volume AV fluctuations are not statistically 
independent, the bulk modulus is related to the matrix 
[C] via V/kt = N^^^ + A^t^iCn + ^i^tCit + ^l^ii- 
The positivity of the volumetric specific heat Cv , chemi- 
cal susceptibility matrix [^] and bulk modulus B = 1/kt 
are the necessary and sufficient conditions for thermody- 
namic stability, which must be satisfied irrespective of 
approximations used at the microscopic level. 



1 and US-2 is made via thermodynamic prop- 
erties such as the isothermal compressibility 
KT = {V/N^) {dN+/di-i+)rpy , which is shown in 
Fig. [5] versus l/{kF+as) for the values of the Zee- 
man field hz/eF+ = 0,0.5,1.0,2.0 and spin-orbit 
couphng v/vp^ = 0.28. Notice the negative re- 
gions of kt indicating that the uniform superfluid 
phases are unstable, and its discontinuities at phase 
boundaries. The normal state compressibility kt or 
R.T — {N^)V~^ Kt = {dN-^/dfi+)rp y can be obtained 
analytically for arbitrary Zeeman and spin-orbit 
parameter v in the BCS limit where l/kp^ag 



-oo as 



37V4 
47" 



E 



where the auxiliary function Aj is 



(11) 



A. 



+ V) ± J {fL+ + V) 



hi 



and its derivative is 




16 
12 

+ 

8 
c 

4 
0, 



1 1 


/ 




/ 50| l-n 1 

-50' " ' 

0.5 1 



-10 12 




12 3 4 

1/kF+as 



dA± 
djj.+ 



l±vlJ{lL++vY-{p\-hl) 



/(2A±) 



with = [Lj^jepj^, hz = hz/tp+, and v = v/[2ep^). 
Notice that, as ft,^ — )• and v 0, A± v^mT and 
Kt — >■ (3A^+)/(2ei?+) is reduced to the standard result, 
since p.-^- — t- 1. In addition, kt or kt can be obtained 
analytically in the BEC limit where l/kp^ag — > +00. 
When /i, and v are zero, then 



Kt = 



3N+ TT 

2eF+ kp^Qs 



can also be written in terms of bosonic properties 



1 /dN. 



1 / rriB 



V \ dfj.+ J j,y TT V as 



(12) 



(13) 



FIG. 5: Isothermal compressibility kt = {N\)V ^kt ~ 
{dN+/dfi+)rp Y in units of 3A^+/(4eF+) versus interaction 
l/{kF+as) at spin-orbit coupling v/vp^ — 0.28 for tiie val- 
ues of the Zeeman field a) h^/eF+ ~ 0; b) hz/ep^ = 0.5; 
c) hz/eF+ = 1.0; and d) hz/eF+ = 2.0. Insets show regions 
where the compressibility is large. 



Further characterization of phases US-0, US- 



where niB = 2to is the boson mass and as = 2as in 
the boson-boson interaction. In the case where /i^ 7^ 
and V ^ Q, a, similar expression can be derived for 
{dNj^ / dii+)p y but the effective boson mass ms = 
2mf{hz,v), and the effective boson-boson interaction 
as = 2asg{hz,v) are now functions of hz and v. Notice 
that the ratio tub /as in the BEC limit can be directly 
extracted from the behavior of Rp for large \/{kp^as). 



